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■ Abstract. We analyze 9 open clusters with ages in the range 70 Myr to 3.2 Gyr using UCAC2 proper motion data 

and 2MASS photometry, which allows us to reach stellar masses down to w 0.7 Mg. We employ in this work an 
approach in which the background proper motion contribution is statistically subtracted in order to obtain the 
1 cluster's intrinsic proper motion distribution. For each cluster we consider the projected velocity distributions in 

the core and off-core regions separately. In the projected velocity distribution of all sample clusters we find a well- 
defined low-velocity peak, as well as an excess in the number of stars at larger velocities. The low-velocity peak 
is accounted for by the random motion of the single stars, while the high-velocity excess can be attributed to the 
large velocity changes produced by a significant fraction of unresolved binaries in a cluster. We derive kinematic 
parameters of the single-star distribution, in particular the projected velocity dispersion. The relatively large 
velocity dispersions derived in this work may reflect the non-virialized state of the clusters. Based on the relative 
number of high-velocity (binary) and single stars, we inferred for the sample clusters unresolved binary fractions 
in the range 15% < fbin < 54%, for both core and off-core regions. Stars with a projected velocity exceeding the 
' maximum reached by the single-star distribution are identified in 2MASS J x (J — H) colour-magnitude diagrams. 

The asymmetry observed in the distribution of these stars around the main sequence is consistent with models 
of main-sequence widening resulting from unresolved binaries combined with 2MASS photometric uncertainties. 
The present results suggest that care must be taken when applying proper-motion filters to sort out members, 
especially binaries in a star cluster. This paper shows that proper motions turn out to be a useful tool for 
identifying high- velocity stars as unresolved binary cluster members, and as a consequence, map and quantify the 
binary component in colour-magnitude diagrams. 
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1. Introduction the fraction of multiple systems both increases and tends 

' to concentrate in the central regions, thus changing the 

Binaries play a vital role in the dynamics of star clusters initiai spatial distribution of stars (Takahashi & Portegies 

since, as a result of the combined effect of high star density 2 wart 2000) 
and low velocity dispersion, the energy locked up in binary 

binding energy may be comparable to, or even exceed that Observationally, the main effects of a significant frac- 
of the whole cluster (Hut[1996|). Consequently, changes in tion of unresolved binaries in the central parts of a star 
binary properties along the normal stellar evolution will cluster are expected to be: (i) the production of large de- 
affect the dynamical evolution of the cluster as a whole, viations in the cluster's proper motion distribution with 
particularly the old, more massive clusters. In this context, respect to the single-star component, since the resulting 
the identification and quantification of the binary fraction velocities may be much larger for the binary motions, (ii) 
is fundamental for understanding the dynamical evolution the underestimate of the number of low-mass stars with 
of star clusters. respect to higher-mass stars, and (Hi) the widening of 
Open clusters are known to contain significant frac- the main-sequence (MS) towards higher magnitudes in a 
tions of the stellar content in the form of multiple sys- colour-magnitude diagram (CMD). Point (i) is the fun- 
tems, binaries in particular (e.g. Montgomery, Marschall damental issue addressed for the first time in the present 
& Janes |1993 von Hippel & Saraiedini ll998|l . And, as a paper. As a consequence of (ii), the observed central lumi- 
consequence of the dynamical evolution in open clusters, nosity (or mass) function turns out to be flatter than that 
of the actual, resolved stellar distribution. In addition, the 
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binary fraction and MS mass range (e.g. Kroupa, Tout & 
Gilmore }199 II and references therein). 

Rather than consider the proper motion components 
of individual stars, which might exclude unresolved bina- 
ries as non-member stars, we apply in this work a sta- 
tistical analysis. Our approach is to compare the proper 
motion (modulus) distribution of the stars in the cluster 
field with that of the offset field stars. Thus, the presence 
of unresolved binaries should be indicated by significant 
deviations in the cluster's proper motion modulus distri- 
bution. 

In the present work we employ J and H 2MASS 1 data, 
which provide photometry for stars in a large spatial area 
in the direction of a cluster, and proper motion data from 
UCAC2 2 in the same field. 

This paper is organized as follows. In Sect. we provide 
general data on the target clusters. 2MASS photometry, 
CMDs and structural parameters for the clusters are pre- 
sented in Sect.|3| The proper motion analysis is presented 
in Scct.^J Finally, concluding remarks are given in Sect. [5] 



2. The target clusters 

The open clusters selected for the present analysis are 
M26 (NGC6694), NGC2287, M48 (NGC2548), M93 
(NGC2447), NGC5822, NGC6208, NGC3680, IC4651 
and M 67 (NGC 2682) . These clusters span a wide range 
in age, which is important, since the binary fraction in a 
cluster is expected to increase as a function of dynamical 
evolution (age). General data for the above open clus- 
ters are given in Table The data in cols. 2-6 have 
been obtained from the WEBDA open cluster database 
(Mermilliod JT996) 3 . In cols. 7-11 we list parameters de- 
rived from the present work using 2MASS data (Sect.[3Jl. 

In previous studies binaries have already been detected 
in some of the above clusters. Stetson 1)1981(1 indicated 
that NGC 5822 contains a high percentage of binary stars 
and several blue stragglers. Levato & Malaroda l|1979[) 
found that most of the B8-A0 stars in NGC 2287 are 
spectroscopic binaries. Nordstrom, Andersen & Andersen 
( 1997) identified several single and binary cluster members 
in NGC 3680. With respect to IC 4651, Meibom, Andersen 
& Nordstrom PUU2|) found that 37% of the 19 red giants 
are spectroscopic binaries with periods of up to 5 000 days, 
while 52% of the 67 MS stars are binaries with periods less 
than 1000 days. Montgomery, Marschall & Janes (1993) 
found that at least 38% of the stars in M 67 are binaries. 
Mermilliod & Mayor (1982J) have found a percentage of 
25-33% of binaries in M 93. 



1 The Two Micron All Sky Survey, All Sky 
data release (Skrutskie et al. HWj7fl . available at 
http:/ / www. ipac. caltech. edu/2mass /releases/ allsky/ 

* The Second U!s! Naval Observatory CCD 
Astrograph Catalog (Zacharias et al. f2D03|, available at 
http://vizier.u-strasbg.fr/viz-bin/ VizieR?-source=UCAC2 

^ http://obswww.wnige.ch/webda 



3. The 2MASS photometry 

2MASS has proved to be a useful tool to explore vari- 
ous aspects related to open clusters, e.g. Bica, Bonatto, & 
Dutra (2003, 2004}. We derive parameters for the present 
cluster sample from J and H 2MASS photometry. The 
VizieR tool 4 was used to extract stars in a circular area 
centered on the coordinates given in Tabled In order to 
maximize the statistical significance and representativity 
of background star counts, we used stars extracted in an 
external annulus as offset field. 

3.1. Colour-Magnitude Diagrams 

In Fig. ^ we present the J x (J — H) CMDs for the open 
clusters dealt with in this paper. To maximize cluster 
membership probability, the stars included in the CMDs 
correspond to a region with radius 2 x R cor e (Sect. 13. 2|) . 
Cluster parameters were derived by fitting solar metal- 
licity Padova isochrones (Girardi et al. I2002|l computed 
with the 2MASS J and H filters 5 , to the observed CMDs. 
The 2MASS transmission filters produced isochrones very 
similar to the Johnson ones, with differences of at most 
0.01 in (J - H) (Bonatto, Bica & Girardi l2lTO|l . see also 
Grocholski & Sarajedini (2003) for an analysis of combined 
optical/near- infrared CMDs. For reddening and absorp- 
tion transformations we use Ry = 3.2, and the relations 
Aj = 0.276x Ay and E(J — H) = 0.33 x E(B - V), follow- 
ing Dutra, Santiago & Bica ( 2002) and references therein. 
The resulting age, colour excess and distance to the Sun, 
are given respectively in cols. 7, 8 and 9 of Tabled Except 
for M 26 and NGC 6208, no amount of reddening was de- 
tected for the clusters. 

For most clusters, the optical value for E(B — V) in 
TableQlis larger than that inferred from the near-infrared. 
However, we note that E(B — V) = 3 x E(J — H), so these 
near-infrared bands are in fact less sensitive to redden- 
ing determinations. In previous studies we already noticed 
the need for less equivalent reddening in the near-infrared 
than what was expected from the optical, e.g. the case 
of M67 (Bonatto & Bica HSSBt - The excesses E(J - H) 
are derived from the overall isochrone fit to the CMD se- 
quences. Besides, for several of the present clusters, deep 
CMDs are studied. As can be seen in the CMDs of Fig.^ 
there is not much freedom for variations in E(J — H) mea- 
surements. Basically, all values of E(B — V) smaller than 
0.06 are null in the near-infrared. In the case of IC4651 
(Fig. HJ), the WEBDA value E(B - V) = 0.12 would pro- 
hibitively shift the isochrone too far into the red. 

3.2. Cluster structure 

Structural parameters for the clusters have been derived 
by means of the radial distribution of stars, defined as the 

4 http://vizier.u-strasbg.fr/viz-bin/VizieR?-source= 11/246 

5 available to download at 
http:/ /pleiadi.pd. astro, it/isocjphotsys. 01 /isocjphotsys. 01. html 
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243 


0.03 


0.69 


160 ± 10 


0.00 


0.76 ± 0.04 
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8.6 ±0.4 


M48 


08:13:43 


-05:45:00 


360 


0.03 


0.77 


360 ± 40 


0.00 


0.76 ± 0.03 


0.9 ±0.2 


8.9 ±0.4 


M93 


07:44:30 


-23:51:24 


387 


0.05 


1.04 


400 ± 50 


0.00 


1.05 ±0.04 


0.6 ±0.1 


5.5 ±0.6 


NGC 5822 


15:04:21 


-54:23:48 


662 


0.15 


0.92 


1000 ± 100 


0.00 


0.72 ± 0.03 


0.9 ±0.1 


8.5 ±0.4 


NGC 6208 


16:49:28 


-53:43:42 


1170 


0.21 


0.94 


1300 ± 100 


0.21 ±0.03 


1.31 ±0.06 


0.9 ±0.2 


6.5 ±0.8 


NGC 3680 


11:25:38 


-43:14:36 


1200 


0.07 


0.94 


1600 ± 100 


0.00 


1.05 ±0.05 


0.6 ±0.1 


7.0 ± 0.6 


IC4651 


17:24:49 


-49:56:00 


1140 


0.12 


0.89 


1800 ± 200 


0.00 


0.91 ±0.04 


0.9 ±0.1 


5.9 ±0.3 


M67 


08:51:18 


+ 11:48:00 


2560 


0.06 


0.91 


3200 ± 100 


0.00 


0.87 ±0.09 


1.2 ±0.1 


8.8 ±0.9 



Table Notes. The 2MASS colour excess for M26 is E(J - H) = 0.14, and for NGC 6208 E(J - H) = 0.07. Uncertainties in 
columns (7)-(9) are derived from the isochrone fit. 
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Fig. 1. J x (J — H) CMDs for the present open clusters. 
In each panel, the region shown corresponds to a radius 
r = 2 x Rcorc- The solid line shows the best-fit Padova 
isochrone. The derived age is indicated in each panel. 

projected number of stars per area around the cluster cen- 
ter. Before counting stars, we applied colour filters to both 
cluster and corresponding offset field, in order to take ac- 
count of the contamination of the Galaxy. This procedure 
has been previously applied in the analysis of the open 



clusters M67 (Bonatto & Bica l2TMJ) . NGC 188 (Bonatto, 
Bica & Santos Jr. l2TM| and NGC 3680 and NGC 2180 
(Bonatto, Bica & Pavani I2004J) . The radial distribution 
was determined by counting stars inside concentric annuli 
with a step of 1.0' in radius. The background contribution 
level corresponds to the average number of stars included 
in the external annulus (offset field). 

The resulting projected radial density profiles for the 
present clusters are shown in Fig. For a better compari- 
son between clusters, the radii on the abscissas have been 
scaled from arcmin to parsecs using the distances derived 
in Sect. 13.11 The statistical significance of each profile is 
reflected by the la Poisson error bars. 

First-order structural parameters for each cluster are 
derived by fitting the two-parameter King l)1966afl sur- 
face density profile to the background-subtracted radial 
distribution of stars. The two-parameter King model es- 
sentially describes the central region of normal clusters 
rKing !1966bl Trager, King & Diorgovskv lTMty . The fits 
were performed using a non-linear least-squares fit rou- 
tine which uses the errors as weights. The best-fit solu- 
tions are shown in Fig.[2as dashed lines, and the resulting 
core radii are given in col. 10 of Table Considering the 
profile fluctuations with respect to the background level, 
we can define a limiting radius (Rn m ) for each cluster. 
For regions beyond Ri; m the null-contrast between cluster 
and background star density would produce exceedingly 
large Poisson errors and consequently, meaningless results. 
Thus, for practical purposes, we can consider that most of 
the cluster's stars are contained within Rn m . The limiting 
radii are given in col. 11 of Tabled 

4. Proper motion analysis 

Proper motion components in right ascension and declina- 
tion for the cluster stars have been obtained from UCAC2. 
Rather than work separately with proper motion compo- 
nents, we employ in the present work the projected ve- 
locity on the sky v p — \f$ x x cos(<5) 2 ± fij, which corre- 
sponds to the modulus of the proper motion components. 




E. Bica and C. Bonatto: Proper motions in open clusters 

7 



R(pc) 



R(pc) 



R(pc) 



Fig. 2. Projected radial distributions of stars. For each 
cluster, the average background level is shown as a narrow 
shaded area; la Poisson errors are also shown. The dashed 
lines show the two-parameter King model fitted to the 
radial distributions of stars. In each panel, the core radius 
is indicated as a dotted line. For a better comparison, all 
panels are plotted on the same scale. 



Although affected by the systemic motion, asymmetries in 
the proper motion distribution of a cluster may yield infor- 
mation on the internal kinematics. To be consistent with 
the 2MASS analysis, we extracted from UCAC2 proper 
motions for each cluster inside the same area as that used 
to extract 2MASS photometry. Since the UCAC2 cata- 
logue also includes the 2MASS photometry, we verified 
that the correspondence of the two catalogues is nearly 
complete for J < 14.5. This is important since for the in- 
terpretation of proper motion data in the context of deter- 
mining the physical nature of an object one must include 
faint stars as well. 

4.1. M67, a test case 

Because of its relatively high Galactic latitude (b « 32°) 
the field of M 67 (NGC 2682) presents little background 
contamination (Fig. QJ. Thus, M67 is an ideal cluster to 
examine the probable association of proper motion with 
the presence of binary stars. In this sense, M 67 is an inter- 
esting cluster because it is known to contain a significant 
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Fig. 3. Correspondence between the 2MASS photometry 
and UCAC2 data for M67. Both catalogues are nearly 
complete for J < 14.5. The colour filter is shown as dashed 
lines. Left panel: central (R < 10' w 2.5 pc) region; right 
panel: offset field (same area). 

fraction of unresolved binaries (Montgomery, Marschall & 
Janes HMty . In Fig. we show the J x (J-H) CMDs 
of the central region of M 67 (left panel) and the corre- 
sponding (same area) offset field (right panel). Stars with 
2MASS detections are indicated by open circles, while 
those with proper motion determined by UCAC2 are in- 
dicated by dots. For stars with J < 14.5 both catalogues 
are essentially complete. We also show in both panels the 
colour filter used to isolate most M 67 stars from the back- 
ground contamination. Probable M 67 blue stragglers (at 
10 < J < 11 and 0.0 < (J - H) < 0.3) have also been ex- 
cluded from the analysis. The relatively large number of 
stars in M 67 also allows us to study spatial variations in 
the proper motion distribution. Thus, the present analysis 
will be applied to both core and off-core regions. 

Although small, the contamination by Galactic stars 
of the CMD of M 67 must be taken into account to isolate 
the intrinsic proper motion distribution in this cluster. 
We study the properties of the proper motions in M67 
according to the following procedure. First we apply the 
colour filter shown in Fig. [3] to both cluster and offset 
field stars. The filtering process takes into account most 
of the background (mostly Galaxy disk stars) contamina- 
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Fig. 4. Proper motion components in the core (bot- 
tom panels) and off-core (top panels) regions of M67. 
Histograms in the left panels show the number of stars 
per velocity bin for the [i a component; us components are 
in the right panels. The shaded areas correspond to the 
offset field (same area) contribution. The solid line rep- 
resents the intrinsic proper motion distribution in each 
region. 
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Fig. 5. Distribution of projected velocity of stars in M67. 
Histograms in the left panels show the number of stars 
per velocity bin in the core (bottom panel) and off-core 
regions (top panel). The shaded area corresponds to the 
offset field (same area) contribution. Right panels: pro- 
jected velocity distribution. Dotted lines: Gaussian fits to 
the distributions. Note an additional high-velocity peak 
centered at v p sa 130 km s" 1 in the off-core profile of M 67. 



tion, leaving only a residual contamination which will be 
dealt with by the offset field. Next we build histograms 
with the number of stars in velocity bins of 2masyr~ 1 
width, both for the cluster and the offset field. This pro- 
cedure is applied to both core (0.0 < i?(pc) < 1.2) and 
off-core (1.2 < i?(pc) < 8.8) regions. Then, the offset field 
histogram is scaled to match the projected area of both re- 
gions. Finally, the subtraction of the offset field histogram 
from those of the core and off-core regions produces the in- 
trinsic proper motion distribution in M 67. This procedure 
is illustrated in Fig.^for the proper motion components in 
right ascension (/x Q ) and declination (fj-s)- As expected for 
a self-gravitating system, the proper motion distributions 
in the core and off-core regions of M 67 present systematic 
deviations with respect to the offset field distribution. 

Significant deviations in proper motion with respect 
to a Gaussian distribution are seen in both fj, a and 
components. Considering the offset field contribution, the 
excesses may be accounted for by a component of high- 
velocity stars inside M67. In order to characterize these 



deviations in terms of high- velocity stars, we hereafter 
use the projected velocity v p . The units masyr -1 have 
been transformed to kms" 1 using the distances listed in 
col. 9 of Table The above procedure is illustrated in 
Fig. in which the left panels show the projected veloc- 
ity histograms both for the cluster regions and offset-field, 
while the right panels display the intrinsic, offset field- 
subtracted, projected velocity distribution in M 67. The 
statistical significance of this procedure can be estimated 
from the resulting Poisson error bars. 

The projected velocity distributions both in the core 
and off-core regions of M 67 present conspicuous low- 
velocity and high-velocity peaks. The low-velocity peaks 
can be accounted for by the random motion of the single 
stars, superimposed on the cluster's systemic motion. On 
the other hand, the high-velocity peaks, which consider- 
ing the error bars represent real excesses, may be produced 
by unresolved binary systems, in which the presence of a 
secondary changes appreciably the velocity of the primary 
star. Kinematical parameters for the cluster are derived by 
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fitting Gaussian functions N oc exp — 2 J° to the v p dis- 
tributions, where vq is the peak average velocity and a v is 
the projected velocity dispersion. We found for the core re- 
gion a low-velocity peak centered at vq = 41.3±0.7kms~ 1 , 
with a velocity dispersion <r v = 6.0 ± 0.6 kms -1 . The 
high-velocity peak, with smaller amplitude, is centered 
at v = 72.8 ± 2.3kms _1 , with a v = 11.5 ± 1.5 kms -1 . 
A similar situation is observed in the off-core region, 
with the low-velocity peak at vq — 49.8 ± 1.8kms _1 , 
with a v w 7.5 ± 2.3 kms -1 , and the high-velocity one at 
v a ~ 73.0 ±2.2 kms -1 , with a v « 10.1 ± 0.7 kms" 1 . The 
above fits are shown in the right panels of Fig. as dot- 
ted lines. Within the uncertainties, a v for the low- velocity 
peak is the same, both in the core and off-core regions. 
The same conclusion applies to the high- velocity peaks in 
both regions. 

Considering the v p distribution in the core region, stars 
with velocities larger than 70 kms -1 can be considered to 
have deviant proper motions with respect to the single- 
star distribution. We identify and plot these deviant stars 
in the CMD of the core region in Fig. El in which the 
3.2 Gyr isochrone is also plotted for clarity. The high- 
velocity stars are distributed along the MS with an asym- 
metry towards the right side of the isochrone. They oc- 
cupy the locus encompassed by the colour-magnitude lim- 
its expected for binary stars combined with the 2MASS 
photometric uncertainties (thin-solid lines in Fig.[f)|). 

To simulate the effects of unresolved binaries on the 
2MASS CMD, we first generate a population of single 
stars according to a standard (x = 1.35) Salpeter MF 
distribution. To be consistent with the 2MASS CMD of 
M 67 we used single stars with masses lower than the re- 
spective turnoff. According to a pre-defined fraction of 
CMD binaries, we randomly select stars from the original 
distribution and build the specified number of binaries. 
No mass or luminosity bias is introduced in this process 
(Kroupa 2001; Mazeh et al. I2003f) and, consequently the 
probability of a star being selected as a binary member 
depends only on the number frequency of its mass. In 
order to obtain the maximum spread in (J — H) colour 
spanned by the binary distribution we use a binary frac- 
tion of 100%, which means that all objects in the CMD are 
unresolved binary systems. We also include in the simu- 
lation the effects of the 2MASS photometric uncertainties 
in J and H. In 2MASS photometry, the average photo- 
metric uncertainty (e) for a given magnitude can be ob- 
tained following ej = 0.023 + 4.78 x icr 9 e( J / a976 ) and 
e H = 0.025+1.92xl0~ 8 eW a998) (Bonatto, Bica & Santos 
Jr. EOnil). For each J and H magnitude, the photometric 
error is obtained by randomly selecting a value in the in- 
terval (—6j,h < e < +£,/,#), assuming that e follows a 
normal distribution curve. Finally, for each single-mass 
star and binary in the CMD, uncertainties are assigned to 
the J and H magnitudes as described above. The result- 
ing maximum colour spreads are indicated in Fig. The 
high- velocity stars in Fig. do not reach the red border 
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Fig. 6. Stars with projected velocity larger than 70 km s -1 
(filled circles) are superimposed on the CMD of M 67. The 
thin solid lines mark the borders of the CMD region oc- 
cupied by unresolved binaries combined to the 2MASS 
photometric uncertainties. The 3.2 Gyr isochrone is also 
shown for clarity. 

of the simulation. Consequently, the binary fraction must 
be lower than 100%. 

Under the assumption that the low-velocity peak es- 
sentially represents the velocity distribution of the single 
stars, and that the high- velocity excess is caused by unre- 
solved binaries, the binary fraction (fbm) in & cluster can 
be estimated by dividing the number of high- velocity stars 
by the total number of stars. Thus, in the core region of 
M 67, we find fbm = 39 ± 16%, which agrees with the frac- 
tion derived by Montgomery, Marschall & Janes l|1993J) . In 
the off-core region we find fbi n = 45 ± 12% which, within 
the uncertainties, agrees with the binary fraction in the 
core. 

4.2. The remaining sample 

The remaining clusters have been analyzed in a similar 
way as M 67, and the corresponding projected velocity dis- 
tributions for the core and off-core regions are shown in 
Figs. and [SJ Clearly, the low- velocity peak is present 
in both core and off-core regions in all clusters, with 
varying central (systemic) velocity and velocity disper- 
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sion. This is reflected in the parameters of the Gaussian 
fits to the low-velocity peaks, which are summarized in 
Table El We include in Table El the maximum velocity 
reached in the single-star distribution, v cut . The high- 
velocity excess is conspicuously present in most clusters 
as well, the exceptions being NGC3680 and NGC 5822 
(core and off-core regions), and M48 and M93 (core re- 
gion). The large error bars present in the profiles of the 
latter clusters, which reflect a relatively small number of 
stars, make the excess marginal. Gaussian fits have not 
been applied to the high-velocity excess in the remaining 
clusters because of the less well defined profiles. Besides, 
as a consequence of the large error bars the resulting 
binary fractions contain significant uncertainties. Within 
the uncertainties, the average binary fraction in the core, 
< fbin >= 35 ± 1 5%, turns out to be the same as in the 
off-core region, < fbin >= 36 ± 11% (standard deviation 
from the average). The same is valid for the average ve- 
locity dispersions in the core, < cr v >— 8 ± 4kms~ 1 , and 
off-core regions, < a v >= 9 ± 2kms _1 . The high- velocity 
stars (according to col. 10 of Table EJ) that are present in 
the core region of each of the remaining clusters are iden- 
tified and plotted on the corresponding CMDs in Fig. [9] 
For each cluster we simulated the combined effect of bi- 
naries and photometric errors, following the approach ap- 
plied to M67 fSect. l4~TJl . The resulting limiting lines are 
shown in the figure. In most clusters the asymmetry in 
colour presented by the high-velocity stars with respect 
to the corresponding isochrone is similar to that displayed 
by the core of M67 (Fig. |HJ), bearing in mind that M67 
is more populous. The observed asymmetries in the CMD 
distributions support the scenario of unresolved binaries, 
as further indicated by the present proper motion study. 

In the CMD of some clusters, a nearly vertical red 
sequence shows up (e.g. IC4651) which we suspect to be 
solar neighbourhood high-velocity low-mass stars. From 
their (J — H) colours, we infer them to be K and M dwarfs. 
Similar to the case of M 67, most of the high- velocity stars 
do fall within the red border of the binary simulation, 
except in M 26 and NGC 6208 which are much affected by 
red stars. 

According to the virial theorem, the dynamical mass 
for a gravitationally bound, spherically symmetric and 
isotropic cluster is given by M = 3a 2 R/G, where a is the 
velocity dispersion and R is the effective gravitational ra- 
dius. As a first-order approximation, we assume R = R core 
and a — a v (core) . For the range of core radii (Table ^) 
and velocity dispersions (Table EJ found in the present 
cluster sample the implied dynamical masses would be in 
the range (1 — 4) x 10 4 M© . For M 26, the youngest cluster 
in our sample with the largest o~ v , a mass would result 
of ~ 2 x 10 5 Mq. These mass values turn out to be ex- 
ceedingly high for open clusters, since the old (« 7 Gyr) 
and very populous open cluster NGC 188 has an estimated 
mass of « 4 x 10 3 M©, as determined from different meth- 
ods (Bonatto, Bica, & Santos Jr. 120041 and references 
therein). The present clusters are clearly not as populous 
as NGC 188. Consequently, the discrepancy in the mass 
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Fig. 7. Projected velocity distributions in the core (left 
panels) and off-core (right panels) regions. The dashed 
curves represent a fit to the single-star distributions, ac- 
cording to the parameters in Table [3 

determinations certainly reflects the non-virialization of 
the clusters. 

The present sample contains one cluster (NGC 2287) 
in common with the clusters studied in Sagar & Bhatt 
(1989). Their analysis was based on radial velocity and 
proper motion measurements of individual stars, re- 
stricted to the brighter, higher-mass cluster members. In 
contrast, the present 2MASS photometry depth allowed 
us to reach low-mass stars as well, as can be seen in col. 2 
of Table El 

Under the first-order assumption of energy equiparti- 
tion and different mass groups within a cluster, the ve- 
locity dispersion would follow the relation a ~ 1/ \fwi, 
where m is the average mass of each stellar group. Thus, 
IIM 



Sagar & Bhatt's (1989) results should favour lower veloc- 
ity dispersions than the present work. Indeed, the velocity 
dispersions derived in the present work for the core and 
off-core regions of NGC 2287 are larger than those found 
by Sagar & Bhatt (1989) in similar regions, by a factor of 
w 6 in the core and « 3.4 in the off-core region. The sig- 
nificant differences in velocity dispersion in both studies 
may be accounted for by the larger number of stars, ~ 357 

- including stars as faint as J sw 15 and mass m > 0.6 M© 

- in the present work, as compared to rs 65 (m > 1.7 M©) 
stars in Sagar & Bhatt's (1989) sample. The large num- 
ber of stars allowed us to fit the appropriate Gaussian to 
the background-subtracted projected velocity distribution 
of NGC 2287 (Fig. Ej). Thus, for such a relatively young 
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Table 2. Properties of the single-star distribution and binaries. 
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Table Notes. mi ow is the stellar low-mass end reached in this work; vq and a v have been obtained from the Gaussian fit to the 



low-velocity distribution; fbin is the binary fraction; v cut is the maximum velocity reached by the single-star distribution. 
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cluster as NGC 2287, in which the low-mass stars are ex- 
pected to be far from virialized, the velocity dispersion 
derived from a large mass range should be larger than 
that obtained from higher-mass stars. In addition, Sagar 
& Bhatt ( 1982J) derive, for the very young, non- virialized 
open clusters NGC 2669 and NGC 4755, velocity disper- 
sions of the order of lOkrns -1 . Interestingly, the velocity 
dispersion in the core of M26, the youngest cluster in the 
present sample (« 70Myr), turns out to be the largest, 
which agrees with strong non-virialization. 



Fig. 9. Distribution of the high-velocity stars in the core 
region of the remaining clusters. The minimum value for 
a high- velocity star (v cu t) is indicated in each panel. For 
ease of visualization, only the isochrone and the borders 
of the binary±photometric error simulation are shown. 

5. Concluding remarks 

In this paper we have analyzed 9 Galactic open clusters, 
spanning a range of ages from 70 Myr to 3.2 Gyr, in terms 
of UCAC2 proper motion data and 2MASS J and H pho- 
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tometry. Considering the depth of the 2MASS photome- 
try for the stars with UCAC2 data, the stellar low-mass 
end in the present analysis reached down to w 0.7 M . 
Rather than considering individual stars, we applied in 
this work a statistical approach in which the cluster's 
intrinsic proper motion distribution was obtained taking 
into account the background proper motion contribution. 
Our goals were to infer the intrinsic binary fraction in 
the clusters and derive kinematic parameters of the sin- 
gle stars, the line-of-sight velocity dispersion in particu- 
lar. The large velocity dispersions derived in the present 
work may reflect the non-virialization of the clusters, par- 
ticularly with respect to the low-mass stellar component. 
For each cluster we considered separately the projected 

velocity (v p — y/ p? a x cos(<5) 2 + jti|J distributions in the 

core and off-core regions. We found in all clusters v p 
distributions which are characterized by a conspicuous 
low-velocity peak, as well as an excess in the number of 
stars for larger velocities, in varying proportions for dif- 
ferent clusters. The low-velocity peak can be accounted 
for by the random motion of the single stars superim- 
posed on the cluster's systemic motion. The high-velocity 
excess is attributed to the significant velocity changes 
produced by the presence of unresolved binaries in the 
clusters. The binary fraction is inferred by taking into 
account the relative number of high-velocity and single 
stars. We derived a fraction of unresolved binaries in the 
range 15% < fbi n < 54%, for both core and off-core re- 
gions in the sample clusters. We showed that the asym- 
metry observed in the distribution of high- velocity stars in 
2MASS J x (J - H) CMDs is consistent with models of MS 
widening resulting from unresolved binaries associated to 
2MASS photometric uncertainties. 

A possible application of the present method could be 
in a discussion of the low binary fraction found in the 
globular cluster M4 (1-2%), a surprising result as pointed 
out by Richer et al. (2004), when compared to previous 
estimates for globular clusters, which amount to < 10% 
(Richer et al. HMTjl . Richer et al. applied to M4 

a proper-motion filtering using multi-epoch HST images 
over 6 years. Possibly in this method a significant fraction 
of binaries was overlooked because large-proper-motion 
stars were discarded. 

We conclude that proper motions are a useful tool for 
identifying high- velocity stars as unresolved binary cluster 
members, and as a consequence, mapping and quantifying 
the binary component in CMDs. 
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